Band gap in graphene induced by vacuum fluctuations 
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The electrons in undoped graphene behave as massless Dirac fermions. Therefore graphene can 
serve as an unique condensed-matter laboratory for the study of various relativistic effects, including 
quantum electrodynamics (QED) phenomena. Although theoretical models describing electronic 
properties of graphene have been elaborated in details, the QED effects were usually neglected. 
In this paper we demonstrate theoretically that QED can drastically modify electronic properties 
of graphene. We predict the following QED effect — the opening of the band gap in a graphene 
monolayer placed inside a planar microcavity filled with an optically active media. We show that 
this phenomenon occurs due to the vacuum fluctuations of the electromagnetic field and is similar 
to such a well-known phenomenon as a vacuum- induced splitting of atomic levels (the Lamb shift). 
We estimate the characteristic value of the band gap and find that it can sufficiently exceed the 
value of the Lamb shift. 



PACS numbers: 78.67.Wj, 31.30.jf 
I. INTRODUCTION 

Graphene — a monolayer of carbon atoms — pos- 
sesses unusual physical properties that make it attrac- 
tive for various applications Usually treated as a plat- 
form for the novel high-speed electronics!^ graphene is 
of great interest from the point of view of the fundamen- 
tal physics as well. Indeed, the low-energy electron exci- 
tations in graphene are massless Dirac fermions with the 
linear energy spectrum, e(k) = | k | EHU That makes 
graphene a condensed-matter playground for the study 
of various relati vistic quantum phenomena, such as the 
Klein tunnellingSSEH an d the Casimir effect PHU Up to 
now, most of graphene-related studies were focused on 
its unusual transport properties, and quantum electrody- 
namics (QED) effects arising from interaction of electrons 
in graphene with a quantized electromagnetic field were 
ncg lected.E This paper is aimed to fill partially this gap 
in the theory. We show that due to the giant Fermi ve- 
locity of electrons in graphene, vp ss c/300, QED effects 
are pronounced and can lead to qualitative modifications 
of the spectrum of elementary excitations. 

The linear energy spectrum of electrons in graphene 
comes from its specific honeycomb lattice structure which 
makes the band gap between the valence and conductiv- 
ity bands to be exactly zero.^ There is the long-standing 
problem of the opening of the band gap. The appear- 
ance of a controllable band gap is required for various 
electronical and optical applications of graphene! 1 * 15 * 16 ! 
Aside from this, it is interesting from the fundamental 
viewpoint to analyze how massless Dirac fermions can 
acquire a mass. This question is relevant, particularly, 
in the context of the observation of Majorana fermions 
in condensed matter systems.^ Several mechanisms of 
the band gap opening in monolayer graphene have been 
proposed. Among them are breaking of the symme- 
try between two sublattices of the honeycomb lattice of 



graphene) 18 * 19 * the spin-orbit coupling 2 ^ and the many- 
body interactions leading to the excitonic instability! 21 * 22 * 

Recently, one of us put forward the proposal of open- 
ing the band gap by illuminating graphene with a circu- 
larly polarized light. ^ In this case the gap in the spec- 
trum of elementary electron excitations appears due to 
the formation of composite electron-photon states which 
are similar to polaritons in ionic crystals and quantum 
microcavitiesp3H23 ft should be noted that, within the 
framework of QED, the electron-photon interaction can 
be observed even if "real" photons are absent and elec- 
trons interact only with vacuum fluctuations of electro- 
magnetic field due to emitting and reabsorbing virtual 
photonsPSl Therefore, one can expect that the photon- 
induced splitting of valence and conductivity bands in 
graphene^ will take place due to the vacuum fluctuations 
even in the absence of an external field pumping. This 
QED effect is similar to the well-known Lamb shift in the 
atomic physics, i.e., the vacuum-induced splitting of the 
states 2s 1/2 and 2pi/ 2 of a hydrogen atom with the char- 
acteristic splitting energy A w 4 fj,eV. The Lamb shift, 
discovered experimentally by Lamb and RetherforcpS* and 
theoretically explained by Bethe^ more than 60 years 
ago, is extremely important for understanding and verifi- 
cation of basic principles of QED. That is why it attracts 
the undivided attention of the physics community up to 
nowF^ 

Since clockwise and counterclockwise polarized pho- 
tons shift electron levels in graphene in mutually oppo- 
site directions,^ the band gap opening needs breaking 
the symmetry between virtual photons with different cir- 
cular polarizations. This can be achieved by placing a 
graphene monolayer inside a planar cavity filled with an 
optically active material (see Fig. [I]). As it will be shown 
below, in this case the vacuum fluctuations lead to the 
opening of the band gap in graphene even in the ab- 
sence of an external circularly polarized optical pumping. 
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FIG. 1: (color online) Sketch of the system. A graphene 
sample placed inside a planar cavity filled with an optically 
active media. The arrows with signs + and — correspond 
to clockwise and counterclockwise circularly polarized virtual 
photons, respectively, which are emitted and reabsorbed by 
electrons in graphene. 

It should be noted, that the QED mass renormalization 
in an optically active media has been considered,^ but, 
surprisingly, the most interesting case of massless Dirac 
fermions was not analyzed before. 

II. THE MODEL 

Let us consider the problem of interaction between a 
single electron in graphene and a single photon mode 



of a planar microcavity. Generally, electron states in 
graphene near the Fermi energy are described by the 
eight-component wave function which accounts for two 
elementary sublattices of graphene, two electron valleys, 
and two orientations of electron spinP In what follows 
intervalley scattering processes and spin-flip effects will 
be beyond consideration, which reduces the number of 
necessary components of wavefunction to two. 

The single-particle Hamiltonian of electron in graphene 
coupled to the cavity mode reads (see the Appendix A 
for details of the derivation) 

Ti = 'Hfield + + %int > (1) 

where 

Afield = y^^q, ± «q, + Qq,± ( 2 ) 

q,± 

is the photonic part of the Hamiltonian written in the 
basis of circularly polarized states, 

Hk = hvpfT ■ k (3) 

is the electron Hamiltonian near the point where the 
valence and conductivity bands of graphene touch each 
other (the Dirac point) and 
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is the Hamiltonian of electron-photon interaction in the 
cavity. For defmiteness, we assume the graphene sheet to 
be placed in the center of the cavity. In Eqs. ([2])-(|4| the 
subscript indices, ±, correspond to the photon modes 
with clockwise and counterclockwise circular polariza- 
tions, k = e x k x + e y k y and q = e x q x + e y q y denote in- 
plane electron and photon wave vectors, respectively, e Xt y 
are unit vectors directed along the x, y-axis, e is the elec- 
tron charge, eo is the vacuum permittivity, L is the dis- 
tance between two mirrors of the planar cavity (the cavity 
length), S is the area of graphene sample, w±, q are the 
eigenfrequencies of clockwise and counterclockwise circu- 
larly polarized photons, <z±, q and q are photonic anni- 
hilation and creation operators. The Pauli vector opera- 
tor, <r, acts in the space of two orthogonal electron states, 
|±), corresponding to the two sublattices of graphene in 
accordance with the following rules: <t z |±) = ±|±) and 
o^l^F) = |±), where a ± — (a x ± ia y )/2. Thus, it corre- 
sponds to the pseudospin of electron. 



Eigenstates of the electron Hamiltonian ([3]) are given 
by the expression^ 



|k,±> = 




I-)) , (5) 



where #k = arctan(fc a //c 2; ) and the signs ± correspond to 
electron states in the conductivity and valence bands of 
graphene (the upper and lower Dirac cones, respectively). 
The corresponding eigenenergies are E2? k = ifiuf |k|. 
The eigenstates of the photon Hamiltonian ([2| can be 
written as |JV± q), where N± „ are photon occupations 
number for photons with different circular polarizations 
(±) and wave vectors q. Then eigenstates of the full 
electron-photon Hamiltonian ([I]) can be decomposed in 
the basis of the orthogonal electron-photon states 

|k,±,AT +) AT_) = | k ,±) (g |JV + , q > ® |JV_, q /) (6) 

with the energies 

£ k°±,Ar + ,Af_ =±hv F \k\+hw +t(l N +tCl + nw_^N_ t(1 i . (7) 
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In order to find eigenstates and eigenenergies of the full 
Hamiltonian ([I]), we will use the perturbation theory, 
considering the interaction Hamiltonian Q as a pertur- 
bation. To calculate the energy corrections in the lowest 
order of the perturbation, one needs to find the eigen- 
values of the 2x2 matrix US 1 ) having matrix elements 
H^)(k) = (k,s,0,0|"H int |k,s',0,0). However, it is easy 
to see that all matrix elements of this type are zero and 
one needs to use the second order of the perturbation 
theory. Physically, we need to account for the follow- 
ing processes: the electron with a wave vector k emits a 
virtual photon with a momentum q and then reabsorbs 
it. Note, that in such a process the momentum of the 
electron in the initial state, k, should be equal to its mo- 
mentum in the final state, but the value of the index s 
can be changed: the electron can remain in the same 
Dirac cone or move from one Dirac cone to another one. 
The last process becomes efficient around k = point 
where the energies of the Dirac cones are close to each 
other, which can lead to the lifting of the degeneracy 
as we show below. Therefore, to calculate the spectrum 
of the Hamiltonian ([I]), we need to use the perturbation 
theory for degenerate states. 

Let us briefly remind how the second-order correc- 
tions can be accounted for within the framework of per- 
turbation theory for degenerate states (the details can 
be found, e.g., in Ref. [33]). Imagine that we have 
a set of states {m} which are close in energy to each 

other (this means that the energy distances |em^ — e^,| 
between them are comparable or smaller relative to a 
characteristic energy of the perturbation). The pertur- 
bation does not couple any states to and mf directly 



(otherwise the standard first-order perturbation theory 
is applicable), but couples them to a set of the states 
{/} whose energies lie far from energies of the states 

{to} (this means that the energy distances \e^m — ef' | 
are large as compared with the characteristic energy of 
the perturbation). In our case the set {to} consists 
of the two states {|k, +,0,0), |k, — ,0,0)}, and the set 
{/} corresponds to the states |k', ±, N+ n „, AT_ )71 „/) with 



A+ 



+ N-. 



7^ 0. Then energies of the perturbed 



{to} states can be obtained by diagonalization of the ma- 
trix Hamiltonian 



h = h 



n 



(2) 



(8) 



where H^ is the matrix of unperturbed Hamiltonian ([3j) 
written in the subspace of states {to}, and the matrix 
elements of the Hamiltonian "H*- 2 -* can be found as 



H 



(2) 



x (m\U in t\l}(l\H in t\m'\ 



(9) 



where the summation goes over all set of the states {I}. 

— (2) 

In the case we consider, the Hamiltonian for a given 
electron wave vector k can be written as the 2x2 matrix 



H 



(2) 



hv F k + 



-hvFk + H, 



(10) 



where the vacuum-fluctuation corrections H J + , H k 
and 'KZ\~ = ("Hk + ) are given by 



WA,q - V F \k\ 



e 2 V F y—^ f 2 

q WA ,q[(wA, q -«F|k|) 2 -^|k-q|2] ' 



(11) 



k 4ire L 



WA,q + V F \k\ 



^A.q [(^A.q + v F |k|) 2 - v 2 F \k - q| 2 ] ' 



(12) 



7 A,q "F 



4(|k| 2 + |k-q| 2 ) 



47r£oL A~i J [K,q + «F|k|) 2 ~4|k~q|2].[( WA , q -^|k|)2- u 2 ? |k-q|2] ' 



(13) 



and the symbol A = ± corresponds to the two different 
circular polarizations of virtual photons. The physical 
meaning of the terms in the Hamiltonian ( 10 1 is the fol- 
lowing. The matrix element Ti^ + corresponds to the pro- 
cess: an electron in the upper Dirac cone emits a virtual 
photon and then reabsorbs this photon while returning 
to the same cone. The matrix element corresponds 
to the same process for the electron in the lower Dirac 



cone. The off-diagonal matrix elements 



= HZ 



correspond to the processes in which the electron after 
reabsorption of the photon changes the Dirac cone. Di- 
agrammatic representation of these terms is shown in 
Fig. [2] Diagonalization of the Hamiltonian ( 10 ) gives 



the renormalized energy spectrum of the elementary ex- 
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FIG. 2: Diagrammatic representation of the terms entering in 
the Hamiltonian ( ToJ) . The solid lines correspond to the elec- 
trons and the dashed lines correspond to the virtual photons. 
Index A = ± denotes the two different circular polarizations 
of the photons, and the index s — ± denoted the two different 
Dirac cones. Summation should be performed over both the 
indices A and s. 



In the discussion above we restricted our analysis to 
the single-mode approximation, accounting for the cou- 
pling of the electron in graphene with only one photon 
mode. Going beyond this approximation, one needs to 
perform the summation over all modes n in Eqs. (11) 
(13 1. Keeping in mind that photon modes with even 



numbers n correspond to the zero field intensity in the 
center of the cavity and, thus, do not interact with the 
graphene sheet, one gets the following expression for the 
band gap: 



2tt€qL 



EE i 

n=0 A=± 



XB 



PI 



In 



7T 2 (2n + l) 2 



(18) 

where B± = vf /c± , qo ~ l/&o is the cut-off parameter of 
integration in Eq. ( 13 ) with do being the lattice constant 
of graphene. 



citations in graphene, 



III. RESULTS AND DISCUSSION 



e ± (k)=H+ + /2+H^-/2 



±\ m + /2-H^-/2 + hv F \k\ 



(14) 



Taking into account that at k = we have T~Li 



u-. 



the renormalized electron energy ( 14 ) at the Dirac point 
can be written as 



e ± (0)=H+ + ± %t 



(15) 



It follows from Eq. ( 15 ) that the vacuum fluctuations of 



electromagnetic field in the cavity can open the band gap 
between the conductivity and valence bands of graphene 
at the Dirac point, which is 



£ g — 2 



(16) 



It should be stressed that in the absence of an opti- 
cally active media, the eigenfrequencies of clockwise and 
counterclockwise circularly polarized photons are equal, 
w+, q = c<j_, q . According to the equation ( fi"3|), i n this case 
the term Hq is zero and the band gap ([16]) vanishes. 
Therefore for the band gap opening one needs to fill the 
cavity by an optically active media which splits modes of 
virtual photons with different circular polarizations. In 
this case the photonic dispersions read as 



W ±.q 



(17) 



where q z = ttti/L is the quantized z-component of photon 
wave vector in the cavity, n is the number of photon 
mode, c± = c/n± are the speeds of light with clockwise 
and counterclockwise circular polarizations, and n± are 
the refractive indices for clockwise and counterclockwise 
polarized light, which are different in the optically active 
media, n + ^ n_ . 



It is seen that Eq. ( 18 1 contains the summation over 
photon polarizations, A = ±. As expected, the contri- 
butions of clockwise and counterclockwise polarized pho- 
tons in the band gap ( 18 ) have opposite signs and the 
band gap vanishes for B + = B_ and appears only in 
the presence of an optically active media with B + ^= 8- . 
For instance, the cavity can be filled with a magneto- 
gyrotropic media based on ferrite garnets which possess 
the giant difference between the velocities of light with 
different circular polarizations.^ 4 - The effect becomes even 
more pronounced if the cavity is filled with an active me- 
dia with the circular dichroismPS In this case one of the 
two circularly polarized photon modes in the cavity is 



suppressed and its contribution to the band gap ( 18 1 can 
be neglected. As a result, the summation over A in Eq. 



( 18 1 can be omitted, which leads to the drastic increasing 
of the band gap. 

Figure [3ja) shows the dependence of the band gap, 
e g , on the cavity length, L, for the cavity filled by such 
a media with the circular dichroism. In the physically 
relevant region of the cavity lengths (the white area in 
Fig. [3^a,) ) , the band gap calculated in the single-mode 
approximation (n — 1) is of several /zeV, which is compa- 
rable with the Lamb shift p2 The summation over higher- 
modes increases this value by 1-2 orders of magnitude: 
for 100 modes the value of the band gap increases to 50- 
100 /ieV, while the summation over all modes gives value 
of about 200 /ieV. However, the summation over infinite 
number of modes overestimates the band gap. Indeed, 
if the characteristic photon wavelength, 2n/q z ~ 2L/n, 
is comparable to the interatomic distance, the macro- 
scopic description of an optically active media becomes 
irrelevant. Therefore the band gap can be reasonably es- 
timated to be about tens of ^eV, that is one order of 
magnitude bigger than the Lamb shiftP^l 

The energy spectrum of massive Dirac fermions in 
graphene is plotted in Fig. [3jb) . The renormalized dis- 
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FIG. 3: (color online) (a) The band gap in graphene induced 
by vacuum fluctuations, calculated with accounting different 
numbers of cavity modes n; (b) Energy spectrum of free elec- 
trons in graphene (dashed lines) and electrons dressed by vir- 
tual photons (solid lines). The calculation is performed for 
the cavity length L — 300 nm and the number of accounted 
cavity modes is n = 100. 



persion relation can be approximated by the analyti 
expression 

:vW|k|) 2 + (m*4) 2 , (: 



(k) = ±A 



where the effective mass of electron dressed by virtual 
photons is m* — e g /2vp. 

It should be noted that the considered single-electron 
problem can be easily generalized for the realistic situa- 
tion when the valence band is filled by the Fermi sea of 
electrons. In this case the Pauli principle forbids virtual 
transitions into the lower Dirac cone filled with electrons, 



which reduces both the matrix elements ( 11 1 — ( 13 ) and 



the band gap (18 1 by the factor of 1/2. 



IV. CONCLUSIONS 

We predicted the quantum electrodynamical effect in 
graphene placed inside a planar cavity filled by an op- 
tically active media. Due to the vacuum fluctuations of 
electromagnetic field in the cavity, the spectrum of el- 
ementary excitations in graphene undergoes qualitative 
changes. Namely, the valence and conductivity bands of 
graphene are split at the Dirac points. The value of the 
vacuum-induced band gap can be one order of magnitude 
bigger then the famous Lamb shift in hydrogen atom. 
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Appendix A: Derivation of interaction Hamiltonian 

The introduction of the electron-photon interaction in 
graphene can be done by the conventional replacement 
hk — > Kk — eA, where A is the operator of the vector- 
potential of electromagnetic field. Then the full Hamil- 
tonian of the electron-photon system reads 

H = v F &-(hk-eA S )+~ J dV (e E f eE + /i B f /r ^ , 

(Al) 

where E, B are the operators of electric and magnetic 
fields, and e, fi are the tensors of electric and magnetic 
permittivity of the media, respectively. The integration 
in the last term, giving the energy of free electromagnetic 
field, goes over all space where the field is present. In 
the current paper we consider a graphene sheet placed 
in a planar microcavity. In this case it is convenient 
to represent the operators of the fields in terms of the 
eigenmodes of the cavity as 



A ( r ) = ^™,q( r ) . 

A,n,q 
A,n,q 

B(r) = BA,»,q(r) , 

A,n,q 



(A2) 
(A3) 
(A4) 



where n — 1,2,3, ... is the number of field mode in the 
cavity. Using the Coulomb gauge, we can write the field 
operators (A2)-(A4| as 



L A,n,q 



(r) 



2eoWA,n, 



EA,n,q(r) = i 



2c 



B A,™, q (r) 




aA,n,qUA,n,q(r) + 

+a t A,n, q u L l , q (r)), (A5) 

(4,»,q U A,n,q( r ) ~ 
-flA.n.qUA.n.qWj, (A6) 

aA,«,qV X UA, n ,q(r) + 



+«A,„, q V x u Ajn , q (r) , (A7) 



where u^n.q are the cavity eigenmodes. If the cavity is 
filled with an optically active media, the eigenmodes are 
circularly polarized and can be found as 



where 



u±,„, q (z, r )= e ±V^ sin (^) eiq ' r > ( A8 ) 



w ±,n,q — C± 



7TON 2 



(A9) 



G 



are the photon eigenfrequencies. Therefore the Hamilto- 
nian of the interaction between the graphene sheet and 
the electromagnetic field in the cavity can be written as 

"Hint = -ev F a ■ A(r) = evpV2 ^ (e+o-~ + 

A=±,n,q 

+ e_<7+) • A A ,„, q (r). (A10) 



Since the graphene sheet is placed in the center of the cav- 
ity (at z = L/2), it is coupled only with modes (A8) cor- 



responding to odd numb ers n. This means that the sum- 
mation index, n, in Eq. (A10) is odd: n = 1, 3, 5, 7, . . . . 
Then, using the expression ( A5 ) for the vector potential 



operator A.\ )Jliq (r), the interaction Hamiltonian (A10) 
can be rewritten in the form Q. 
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